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Q . Abstract. By using the recent experimental measurements of B ^ 7I7I and B Kn branching ratios, 

• we find that the amplitudes computed at the leading order of the Aqcd / nif, expansion disagree with 

\ the observed BR&, even taking into account the uncertainties of the input parameters. Beyond the 

i leading order, Charming and GIM penguins allow to reconcile the theoretical predictions with the 
data. Because of these large effects, we conclude, however, that it is not possible, with the present 

^77* ! theoretical and experimental accuracy, to determine the CP violation angle y from these decays. 

' We compare our results with those obtained with the parametrization of the chirally enhanced non- 

' perturbative contributions by BBNS. We also predict large asymmetries for several of the particle- 

^ ! antiparticle B/Js, in pai-ticularB/;(B+ BR{Ba ~* K+n^) and BR{Bd n+%-). 

o : 

^ , 

g : INTRODUCTION 

p • The advent of the B factories, Babar [1] and Belle [2], opens new perspectives for 
p . very precise measurements of non-leptonic 5-decays [3] and calls for a significant 
Qj I improvement of the theoretical predictions. In this respect, important progress has been 

■ recently achieved by systematic studies of factorization in B decays [4, 5]. These studies 
>• . confirmed the physical idea [6] that factorization holds in heavy hadron decays, itiq ^ 

^ ! ^QCD, for the leading terms of the AQco/fnQ expansion. They leave open, however, 
^ I the central question of whether i) the leading terms predict with sufficient accuracy 
the relevant 5-meson decay rates or ii) the power-suppressed corrections, which cannot 
be evaluated without some model assumption, are phenomenologically important. This 
problem has been addressed in a series of papers [7]-[13]. In particular, the main 
conclusion of ref. [12] is that non-perturbative KQcoh^b corrections from the leading 
operators of the effective weak Hamiltonian, conventionally called Charming (or GIM) 
penguins, are very important in cases where the factorized amplitudes are either colour 
or Cabibbo suppressed. One of the consequences of this observation is that factorization 
at the leading order in Agcz) / is unable to reproduce the observed B^rni and B — > Kk 
BRs even taking into account the uncertainties of the input parameters. 

The results of ref. [12] seem to be at variance with the conclusions of ref. [13] where 
it is stated that, in the QCD factorization approach, "an acceptable fit to the branching 
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fractions is obtained even if we impose y < 90° as implied by the standard constraints 
on the unitarity triangle" ^. 

Given the contradictory conclusions of the different studies, it is very important to 
clarify the situation by comparing input parameters, methods of analysis and results. In 
order to help the debate on this subject, we discuss in the following: 

• The theoretical framework of factorization and the calculation of the amplitudes at 
the leading order of the AQco/mh expansion; 

• A comparison of the leading terms, including their uncertainties, with the measured 
B ^ nn and B Kn branching ratios; 

• Results for the 5 — > tcti and B — > Kn BRs including charming (and GIM) penguins; 

• A comparison of the results obtained with charming and GIM penguins [12] with 
the model of the AQco/mh corrections adopted in [13]. 

FACTORIZATION 

In this section we recall the basic ingredients necessary to compute the relevant ampli- 
tudes either whithin factorization or including the corrections which arise at higher order 
in the AgcD/mh expansion. 

The physical amplitudes for B — > Kn and B ^ nn can be conveniently written in 
terms of RG-invariant parameters built using the Wick contractions of the effective 
Hamiltonian [16]. In the heavy quark limit, following the approach of ref. [5], it is 
possible to compute these RG invariant parameters using factorization. The formalism 
of ref. [5] has been developed so that it is also possible to include the perturbative 
corrections to order a,, i.e. at the next-to-leading order in perturbation theory ^. 

For the sake of discussion, it is instructive to start from the explicit expression of the 
B^ — > K'^%~ amplitude. In terms of the parameters defined in [16], this amplitude reads 

siiBd^K+n-) = - VusV:t,(Ei{s,u,u;Bd,K+,n-)-I^^{s,u;Ba,K+,n-)~) 

+ Vt,V;i,Pi{s,u;Bd,K+,%-). (1) 

We have 

Ei{s,u,u;Bd,K+,n-) = fl?(is:7t)(!2i)fact + a2(^^)(G2)fact + ^i 

10 

Pi{s,u;Bd,K+,K-) = £a^'(is:7c)(a-)fact + A 

i=3 

10 

Pr{s,u;B,,K+,%-) = Y^ia^m -a«(i^7t))(a-)fact + ^r , (2) 

i=3 



2 Indeed the analysis of ref. [14] gives Y= (54.8 ±6.2)" and ref. [15] quotes 34" < y < 82". 

^ An alternative framework, provided by the approach of ref. [4], will not be discussed here. For a recent 

discussion see also [17]. 



where (!2()fact denotes the factorized matrix elements and aj the parameters introduced 
in [5] . Eqs . ( 1 ) and (2) are exact. Unfortunately, similar equations require the knowledge 
of several non-perturbative parameters, which at present cannot be extracted from the 
data. To be predictive, we will then use our physical intuition to reduce their number. In 
eqs. (1) and (2): 

1. The CKM matrix elements can either be taken from other experimental measure- 
ments or from a fit to the non-leptonic BRs, assuming that factorization is accurate 
enough; 

2. The coefficients a{{MiM2) (e.g. a"{Kn)) are renormalization group invariant, as 
the corresponding factorized matrix elements, and have been computed perturba- 
tively at the NLO in refs. [5, 13]; 

f f 

3. The coefficients and ag, which have also been computed to one-loop order, are 
instead scheme dependent. Their scheme-dependence is cancelled by the hadronic 
matrix elements of the penguin operators Q(, and Qg respectively. Assuming fac- 
torization for the chirally enhanced contributions [5], the latter can be expressed in 
terms of the ratio 

r (u) = — — — ^ _ , (3) 

^ mb{iJ){ms{ij)+mq{ii))' 

which is formally of O {hQco/tnt) but numerically important (an analogous param- 
eter rJ^ifj) can be defined for %% decays). We will discuss the role of these terms 
below. 

4. The leading amplitudes (2i)fact are computed in terms of decay constants and semi- 
leptonic form factors. The form factors can either be taken from theoretical calcu- 
lations [18, 19] or fitted from the experimental 5/?s (the possibility of extracting 
them from the corresponding B ^ K semileptonic form factor at small momentum 
transfer is at present rather remote). 

5. The tilded parameters, namely Pi, Pf™ and Ei, are genuine, non-perturbative 
^QCD/mb corrections which cannot be computed at present. If we neglect Zweig- 
suppressed contributions, by SU (2) symmetry one can show that all the Cabibbo- 
enhanced Agc£)/w/, corrections toB Kn decays can be reabsorbed in Pi. Several 
corrections are contained in Pi : this parameter includes not only the charming pen- 
guin contributions, but also annihilation and penguin contractions of penguin opera- 
tors. It does not include leading emission amplitudes of penguin operators {Q^-Qiq) 
which have been explicitly evaluated using factorization. Had we included these 
terms, this contribution would exactly correspond to the parameter Pi of ref. [16]. 
The parameter Pi (Pi ) encodes automatically not only the effect of the annihilation 
diagrams considered in [20], but all the other contributions of 0{AQCD/tnf,) with 
the same quantum numbers of the charming penguins. In this respect it is the most 
general parameterization of all the perturbative and non-perturbative contributions 
of the operators and {Q2, and Qa), including the worrying higher-twist in- 
frared divergent contribution to annihilation discussed in ref. [13, 21]. The parame- 
ter Pi has the same quantum numbers and physical effects as the original charming 
penguins proposed in [7], although it has a more general meaning. 



6. If one also includes B — > tuti decays we have several other parameters, for example 
pGiM ^j^^ jj^g formalism of ref. [16]. A closer look to P3 shows that this term 
is due either to Zweig suppressed annihilation diagrams (called CPA and DPA in 
ref. [7]) or to annihilation diagrams which are colour suppressed with respect to 
those entering Pi . For this reason in ref. [12] P3 was taken to be zero. Pi is equal to 
the corresponding parameter in K% decays if SU{3) symmetry is assumed. In our 
analysis we have used the same value of Pi for all Kn and nn channels. Pf ™ will 
be discussed later on. 

We are now ready to discuss factorization for the leading terms of the AQCD/mt, 
expansion. Factorization is the theory of non leptonic decays which is obtained in the 
limit mh — > 00. Thus it consists in neglecting all terms of O {AQco/mh) (Ei = 0, Pi = 0, 

' =0, etc.). At lowest order in perturbation theory, called also naive factorization, the 

a^ are simple combinations of the Wilson coefficients and do not depend on the hadron 
wave functions. The inconvenience of naive factorization is that physical amplitudes still 
have a marked dependence on the renormalization scale because, contrary to the Wilson 
coefficients, the factorized matrix elements are scale independent. 

The scale dependence is reduced by working at O (a^), both for the Wilson coefficients 
and the matrix elements. In refs. [5] it has been shown that, at this order, all the 
dangerous infrared divergences can be reabsorbed in the definition of the hadronic wave 
functions. For the leading terms in the Agczj/w^ expansion there are strong arguments 
to support the idea that this will remain true at all orders in ttv, see also ref. [26]. Thus, 
in the limit nih °°, it is likely that factorization is preserved by strong interactions. 
At 0{as) or higher, the coefficients a{ depend on the specific detail of the hadron 
wave-fuctions. For this reason, the uncertainties relative to the wave functions, as the 
residual renormalization scale dependence, must be taken into account in the evaluation 
of the uncertainties for the theoretical predictions. The approximation in which we 
neglect all the Agco/w^, corrections, but include the perturbative corrections to the 
leading contribution, is called QCD factorization or simply factorization. Factorization 
implies an important consequence: predictions of non-leptonic decay rates are model 
independent to the extent that the few relevant hadronic parameters, namely the kaon 
and pion decay constants, fK,%, the semileptonic form factors, fK,n{0), and the hadronic 
wave functions are known. 

In cases like B Kn decays, where the factorized amplitudes are Cabibbo suppressed, 
the corrections of 0{AQCD/tni,), which unfortunately are model dependent, become 
important. At lowest order in a^, the chirally enhanced terms proportional to (r^) are 
computable by assuming that factorization can be applied beyond the leading order. A 
substantial difficulty arises, however, at O {(XshQCD/mi,). Although the chirally enhanced 
corrections from 06,8 are infrared finite, other contributions of the same order from 
different operators are infrared divergent, signaling that they belong to the class of 
the non-perturbative contributions which appear beyond factorization. These cannot 
be predicted using the same hadronic quantities of the factorized amplitudes. For this 
reason, any phenomenological analysis which aims at including in a coeherent way the 
terms of O {ds^QCD I nih) is forced to introduce extra model-dependent non-perturbative 
parameters besides fK,n{0) and the hadronic wave functions. This implies that, at 



0{KQCD/mi,), model dependence is unavoidable (even in the subsector of the chirally 
enhanced contributions) and it is present in both the analyses of ref. [12] and ref. [13], 
which we will compare below. 

Different model-dependent assumptions were made in the two approaches: 

1. In ref. [12], E\ and P3 were neglected and SU (2) symmetry was assumed for Pi 
(SU (3) when it was used for B ^ %% decays). The same approximations were 
made for P^™. The complex parameter Pi (P^^) was then fitted to reproduce the 
experimental BRs. 

2. In ref. [13] the effects of the chirally enhanced AQco/mt, corrections were either 
computed perturbatively or encoded in two complex phenomenological parameters 
called Xh and Xa. An uncertainty of 100 % to the "default" values (e.g. Xh = 2.4) 
was assigned to these parameters in order to determine allowed bands for the 
predicted BRs. The bands include all other sources of uncertainties. 

We now compare the leading amplitudes to the experimental results in order to test 
factorization. Before using predictions based on factorization to test the Standard Model 
and look for signals of new physics, it is crucial to check how large are the errors induced 
by our ignorance of the O (AQCD/mb) corrections which we are unable to compute. Our 
position, indeed, is that we have more confidence in the SM rather than in factorization. 
In order to test factorization, we ought to use all the information that we have from other 
measurements. Thus, for example, whereas the size of the error on \Vcb\ and \Vub\ can 
be debated, there is no question that these experimental inputs must be included in any 
analysis that aims at testing (or using) factorization. We also stress that the value of | 
is not expected to be affected by the presence of new physics beyond the SM. 

The CP parameter y does in general change if there is physics beyond the SM. It 
remains an interesting exercise, however, to verify whether, by taking the value of y from 
the Unitarity Triangle Analysis (UTA) in the SM, the predicted BRs are in agreement 
with the data. If, by using y from the UTA, one is unable to reproduce the experimental 
B and B — > ^71 BRs, this implies that "either there is new physics or AQco/tnb 

corrections are important" [12]. 

In our analysis we have used the likelihood method which has been described in all 
details in ref. [14]. Without entering in the "ideological" controversy about frequentistic 
and bayesian methods, we only note here that in [14] it has been shown that, at 95 % 
C.L., the Bayesian analysis give the same results as the frequentistic Babar Scanning 
method (and its variations) when the same inputs are used. Thus we will present our 
contour plots, corresponding to fig. 17 of ref. [13], both with factorization plus chirally 
enhanced contributions and with the non-factorizable charming (and GIM) penguin 
corrections. Besides this, we will also give tables with the relevant BRs, both in the 
factorization approximation with chirally enhanced terms and with the charming (and 
GIM) penguin corrections included. 

We end this section with some remarks. In our approach, we have first checked 
that, within factorization and the SM, it is impossible to fit the experimental BRs. 
The AQco/mb terms, that we are then forced to include in order to reproduce the 
experimental results, are non-perturbative quantities, infrared divergent in perturbation 
theory, on which we do not have any knowledge a priori. For this reason we decided to 



fit them on the data. The experimental numbers are nicely reproduced and the corrections 
to factorization are well consistent with the expected size (i.e. Pi is of O (Agco/'w^,) with 
respect to the leading contributions). 

In ref. [13] the subleading power corrections are varied in predefined intervals and 
the change in the predicted BRs is interpreted as uncertainty on the factorized ampli- 
tudes. In our opinion the uncertanties on factorization are only those coming from the 
CKM matrix elements or the form factors, etc. The AQCD/m/, terms instead are really 
contributions beyond factorization: if they are necessary to reproduce the data then it is 
not possible to make model independent predictions. This would remain true even if we 
knew without any uncertainty the hadronic parameters entering at the lowest order of 
the AQco/mh expansion ^. 

Thus we are in the Bermuda triangle: i) without the O (r^^) and O {oisfx) terms, that is 
within QCD factorization, we cannot reproduce the data; ii) the inclusion of the com- 
putable subset of C'(avr^) terms only is inconsistent since there is no reason to exclude 
the other non-perturbative non-computable contributions of the same order. In any case, 
we will show that, by using all the available experimental information, also this case is 
very difficult to reconcile with the data; iii) the complete set of O (cLsr^) corrections leads 
us beyond factorization and the results are model dependent. Indeed there is no proof 
that the one-loop finite chirally enhanced terms remain infrared finite at higher orders in 
Us. Moreover, if the corrections of o{asr'^ ~ a^Agco/w;,) are phenomenologicaly im- 
portant, it is difficult to understand why, at the same level of numerical accuracy, other 
non-chirally enhanced non-perturbative Agco/w^fc terms should not also be taken into 
account. 

The sad conclusion is that the very nice theory of QCD factorization developed in [5] 
is insufficient to fit the data because power corrections, which are model dependent, are 
important in 5 ^ Kn and B ^ nil decays. Finally, model dependence does not implies 
that we are unable to make any prediction. If the assumptions made in our approach are 
reasonable ^, by fitting Pi to the data and with the increasing experimental precision we 
may hope to extract also the value of y or to constrain sin 2a. 



RESULTS 

In this section we present our analysis and a detailed comparison with the results of 
ref. [13]. In order to obtain our results we used the likelihood method as in [12]. The 
input parameters, given in table 1, are also the same but for a few differences: 

• We added the O (a.v) corrections to the coefficients of the penguin and electropen- 
guin operators computed in [13], which appeared after the completion of [12]. In 
this respect the criticisms of ref. [27] do not apply to the present analysis. We will 
see below that, even including these new ingredients, the main physics conclusions 
of ref. [12] are confirmed. 



Note that some BRs are dominated by the AQCD/mi, corrections. 
For example that we may neglect Ei. 



• As for the matrix elements of (26,8 > we include the "computable" factorizable 

chirally enhanced terms in the definition of Pi, which in this way contains all 
the possible Cabibbo enhanced, model-dependent corrections of o{AQCD/tn},). In 

practice this corresponds to fit Pi as in ref. [12] with r^^ ' =0. There is no substantial 
difference between the old choice and the new one, since the contribution of the 
chirally enhanced terms and of Pi have exactly the same quantum numbers. 

• Although this is a minor source of uncertainty we also allow a variation of the 
renormalization scale and of the parameters of the hadronic wave functions. 

We have also verified that by using the inputs of ref. [13] we obtain essentially the same 
results, and hence arrive to the same physics conclusions. 
In this section we present: 

1. A brief discussion of the results obtained in QCD factorization, namely with all the 
terms of O (Agco/w^) set to zero. 

2. The results obtained in our approach by fitting Pi and using y as determined from 
the Unitary Triangle Analysis in ref. [14] 

Y= (54.8 ±6.2)". (4) 

On the basis of this study we predict the B-B asymmetries of the BRs, such as 

^ , BR(B^, ^ %+%-)- BR(B^, ^ %+%-) 

Ji(Ba^Ti+K )^ — ^-4 74 r; (5) 

^ BR{B^^^n+Ti-)+BR(B^^^n+n-) 

3. The results obtained by letting y free and a comparison of our results with those of 
ref. [13]. 



Results in QCD factorization 

In this subsection we compare the model independent results obtained with QCD 
factorization, namely including only the terms which survive when mi, oo, with the 
experimental data. In this case there is still a residual model dependence due to our 
ignorance of the semileptonic form factors fK,n{^) and, at order a^, to the ignorance 
of the hadron wave functions. We vary the semileptonic form factors with flat p.d.f. in 
the intervals given in table 1 and the parameters of the hadron wave functions in the 
intervals given in table 2 of ref. [13]. We also vary the renormalization scale between 
m/,/2 and Imi,, see table 1 of [13]. The data show a generalized disagreement with the 
QCD factorization predictions. In particular the allowed region in the p-T] plane and 
the value of y do not have any overlap with the corresponding ones from the unitarity 
triangle analysis [14]. This remains true even if we double the uncertainty on \Vub\- We 
conclude that QCD factorization cannot be reconciled with data. 



TABLE 1. Input values used in the numerical analysis. The form factors 
are taken from refs. [18, 19], the CKM parameters from ref. [14] and 
the BRs correspond to our average of CLEO, BaBar and Belle results 
[23, 24, 25]. All the BRs are given in units of 10"''. 



A(o) 


0.27 ±0.08 


/if(o)/A(o) 


1.2±0.1 


p 


0.224 ±0.038 


11 


0.317 ±0.040 


BR{Bd K^n'^) 


10.3±2.6 


BR{B+ K+n'^) 


12.0±1.7 


BR{B+ ^K°n+) 


17.4±2.6 


BR{Bd ^ K+K-) 


17.3 ±1.6 


BR{B4^n+n-) 


4.4 ±0.9 


BR{B+ %+%^) 


5.3 ±1.7 



TABLE 2. BRs with Charming ad GIM penguins, with QCD factorization and the infrared- 
finite chirally enhanced corrections only and with the BBNS model which includes Xa.h- All the 
BRs are given in units of 10"*^. 



BR 


Charming 


chirally 


BBNS 


BR 


Charming 


chirally 


BBNS 




+ GIM 


enhanced 






+ GIM 


enhanced 






8.6±0.9 


3.6±1.5 


4.1±1.8 




9.8±1.0 


5.7±2.3 


6.2±2.6 




18.7±1.6 


10.2±4.2 


10.9±4.8 


K+K- 


17.9±1.4 


8.2±3.4 


9.2±4.0 


Jt+7t" 


4.9±0.8 


9.2±3.8 


9.2±3.8 


Tt+Tt" 


3.5±0.8 


5.7±2.2 


6.5±2.5 




0.6 ±0.2 


0.2±0.1 


0.4 ±0.3 











Factorization with Charming and GIM penguins 

We now discuss the effects of charming and GIM penguins, parameterized by Pi and 
P^^. Pi is a complex ampUtude that we fit on the B — > Kn BRs. In order to have a 
reference scale for its size, we introduce a suitable "Bag" parameter, Bi, by writing 

Pl = ^fnW)glBu (6) 

where Gf is the Fermi constant, is a Clebsh-Gordan parameter depending on the final 
Kn (7i7i) channel and B\ = \B\ \ exp(/(j)). Note that B\ differs from the parameter defined 
in ref. [12] because it now includes all the chirally enhanced Agco/w^ corrections, part 
of which were previously explicitly calculated using factorization. In a similar way we 
introduce Bf^. 

We fit all the BRs given in table 1 with GIM and charming penguins included and 
taking the value of y determined from the UTA, see eq. (4). We find 

= 0.13±0.02, (^ = (188±82)^ 

= 0.17 ±0.08, <^^™= (181±59)^, (7) 

where the notation is self-explaining. Note that the size of the charming and GIIVI 
penguin effects is of the expected magnitude. Note that \Bf^\ is very poorly determined. 
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FIGURE 1. p-r| contour plots obtained with factorization and infrared-finite chirally enhanced terms 
using |Vj,f,| (up-left) or letting p and r| free (up-right). We also show a comparison of the p.d.f. of y with 
the one from the UTA analysis of ref. [14] in the two cases. 



The reason is that only BR{B^ — > tt+tt ) in the fit is sensitive to GIM penguins. Thus, in 
practice, we are trying to fit two parameters, namely and to a single BR. 

The results for the BRs can be found in table 2 with the label "Charming-i-GIM". They 
show that the extra charming and GIM parameters radically improve the agreement 
for the measured B Kn and B ^ nn BRs. We do not claim, however, to be able to 
predict BR{Bj tc+tt^) since many effects of the same order besides charming and 
GIM contributions, which in this case are not Cabibbo enhanced, were ignored: our 
results instead show that accurate predictions for Bj rni decays can only be obtained 
by controlling quantitatively the {J^QCD/^b) corrections, which is presently far beyond 
the theoretical reach. To give a complete information, and for comparison with ref. [13] 
we also fit the data by letting y free. In this case we obtain y = (89 ± Al)" . At present, 
the precision of the data and the number of free parameters does not allow a useful 
determination of y. 

The large absolute values of (|), and the sizable effects that penguins have on the 
BRs, stimulated us to consider whether we could find observable particle-antiparticle 
asymmetries as the one defined in eq. (5). We find large effects in BR{B^ K^n^), 



TABLE 3. Absolute values of the rate CP asymmetries 

for B — ^ Kk and B ~^KK decays. 





Charming + GIM 


1^1 


Charming + GIM 




0.05 ±0.03 




0.16 ±0.08 




0.02 ±0.02 


K+n- 


0.15 ±0.07 


7t+7t~ 


0.44 ±0.21 




0.61 ±0.29 



BR{Bd K^n^) and BR{Bd tt+tt^). As discussed before, for BR{Bd n'^n^) our 
predictions suffer from very large uncertainties due to contributions which cannot be 
fixed theoretically. For this reason, the values of the asymmetry reported in table 3 are 
only an indication that a large asymmetry could be observed also in this channel. There is 
a sign ambiguity in ;l ~ sinysincj). This ambiguity can be solved only by an experimental 
measurement or, but this is extremely remote, by a theoretical calculation of the relevant 
amplitudes. For each channel, we give the absolute value of the asymmetry in table 3. 
Note that within factorization all asymmetries would be unobservably small, since the 
strong phase is a perturbative effect of 0(av) [5]. The possibility of observing large 
asymmetries in these decays opens new perspectives. These points will be the subject of 
a future study. 



Comparison with BBNS 

In this section we make a critical comparison with the latest analysis of BBNS [13]. 
We recall that we added all the perturbative corrections and allowed the variations of the 
non-perturbative parameters which were implemented by BBNS. 

We start by discussing the role of chirally enhanced terms that do not suffer from 
infrared divergences at 0{as). This is an istructive case since, if not for other infrared- 
divergent contributions of the same order, the inclusion of only these terms, although not 
justified, would still allow model independent predictions, in the sense discussed before. 
Thus the question is whether these chirally enhanced terms alone, part of which have the 
same effect as Pi, can describe the data. Thus we have repeated the analysis in the UTA 
case with only factorization and chirally enhanced infrared-finite terms with the results 
for the BRs given in table 2 with the label "chirally enhanced". The combined probability 
that all the predicted values for the Kn channels are within two as from the experimental 
numbers is 6%. If we relax the constraint on y we obtain y= (127^=20)" with a 
probability of 0.3%that y < 80'' and of 2% that y < 90^. The corresponding p-T] contour 
plot is given in fig. 1 (upper-left). We conclude that this model (model in the sense of 
including only a chosen subset of the chirally enhanced terms) is strongly disfavoured 
by the data. The reader may be surprised of the difference between fig. 1 and fig. 17 
of ref. [13] (where the non-perturbative parameters Xa,h were however included). The 
difference is explained by the fact that we used the experimental measurements of \Vub\- 
If, following ref. [13, 27], we let both p and r\ free, the p-r| contour plot changes and 
becomes that shown in fig. 1 (upper-right), which is very similar to fig. 17 of ref. [13]. 
For completeness, we also give in fig. 1 the p.d.f. of y in the two cases, together with the 




FIGURE 2. Contour plot in the y- \Vub\ plane, see text. 



p.d.f. obtained with the UTA analysis of ref. [14]. When the experimental information 
on \Vub\ is used the two p.d.f. have no overlap (lower-left). In the other case, it is still 
possible to find values of y compatible with the UTA analysis (lower-right), at the price 
of a rather low value of | Vub \ ■ The situation is well illustrated by fig. 2, where the contour 
plot in the Y-|V„z,| plane of the joint p.d.f. from non-leptonic decays is compared with 
the UTA range for y and the allowed interval for the measured \Vub\, at lo. This figure 
demonstrates that, given the strong correlation between y and |V;,/,|, it is crucial to take 
into account the experimental knowledge of \Vub\ - An important theoretical remark is in 
order at this point. The coefficient is enhanced by the O (a^) corrections and may play 
the same role of charming penguins. It is very scaring that its actual value is strongly 
affected by the contribution of the chromomagnetic operator computed at tree level [13]. 
It is very hard to believe that this contribution which, besides all general considerations, 
is also non-local, will remain infrared-finite in higher orders and can be really evaluated 
in this way. 

We now discuss the effect of Xa^h on the final results. The addition of these param- 
eters, in the range proposed in ref. [13], leaves the situation substantially unaltered as 
shown in fig. 3. This holds true also for the BR%, which are given in table 2 with the 
label BBNS. The reason is that on the one hand Xa and Xh do not have the same quan- 
tum numbers, and hence effects, of charming penguins, on the other the range chosen 
a priori, on the basis on one-loop perturbation theory, is not large enough to improve 
the agreement with the measured BR% (it essentially increases the uncertainty on the 
predictions). Of course by choosing a low value of the strange quark mass and of \Vub\ 
(without using the experimental information coming from its measurements), a large 
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FIGURE 3. p-r| contour plots obtained with the BBNS model using \Vuh\ (up-left) or letting p and 
r| free (up-right). We also show a comparison of the p.d.f. of y with the one from the UTA analysis of 
ref. [14] in the two cases. 



value of /a:(0) and a small value for /7t(0) etc. it is still possible to find some point in 
parameter space where the %^ is good. That this can be used to fit y and test the SM is 
hard to believe though. 

The situation would be different if Xa^h are let free to vary and fitted to the data. We 
have done this exercise and found that the preferred value of is much larger than the 
values allowed in the interval chosen in [13], whereas pn is not determined by the fit. 
We conclude that since Xa and X^ are infrared divergent quantities, the value of which 
cannot be predicted, and since without the inclusion of non-perturbative contributions 
of o(Agc/)/m^,) is not possible to reproduce the experimental data, we are bound to 
use model dependent assumptions in the analysis of non-leptonic B — > Kn and 5 — > tttt 
decays. 



CONCLUSION 



We have analyzed the predictions of QCD factorization for B^rni and B Kn decays. 
Even taking into account the uncertainties of the input parameters, we find that QCD 
factorization is unable to reproduce the observed BRs. The introduction of charming and 
GIM penguins [7] allows to reconcile the theoretical predictions with the data. Instead 
of varying the non-perturbative phenomenological parameters in preassigned ranges, 
we prefer to try to fit them on the data. With the present theoretical and experimental 
accuracy, we find that it is still not possible to determine the CP violation angle y. The 
situation is expected to improve in the near future with more accurate experimental 
measurements of the relevant BRs. 

Contrary to factorization, we predict large asymmetries for several of the particle- 
antiparticle BRs, in particular BR{B+ BR{Bd K^n') and, possibly, 

BR{Bd Ji+Tt^). This opens new perspectives for the study of CP violation in B 
systems. 
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